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ABSTRACT 


The Binary Decision Diagram (BDD) is a very useful representation in the 
design and verification of switching functions. This is due to to its compactness, 
where size is measured by the number of nodes. In the implementation of logic 
circuits, connection of sub-functions is by means of pass transistors. The delay time 
for the interconnections is often larger than the delay of the decision logic. For that 
reason, crossings are often more expensive than logic. Planar Binary Decision 
Diagrams are therefore desirable in implementing logic circuits. This paper presents 
a method for finding a planar Ordered Binary Decision Diagram (OBDD) for 
threshold functions. The program that implements the algorithm is written in 
Borland C++. A special case of Fibonacci threshold function having up to 9 
variables is analyzed. It is shown that Fibonacci functions having up to 9 variables 
have planar OBDD. With this program, the characteristics of other threshold 


functions are developed. 
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Ls INTRODUCTION 

Boolean algebra is used widely in computer science and 
digital system design. Many problems in digital logic design 
and testing, artificial intelligence and combinatorics can be 
expressed in a sequence of operations on Boolean functions. 
However, the classical representation and manipulation of 
Boolean functions have many shortcomings. A variety of methods 
have been developed for representing and manipulating Boolean 
Functions. Those based on classical representations such as 
truth tables, Karnaugh maps, or canonical sum of products form 


are quite impractical--every function of n arguments has a 


representation of size 2” [21 < A more efficient 
representation is the Binary Decision Diagram (BDD). It has 
several advantages. Firstly, most commonly = encountered 


functions have a reasonable representation. For instance, all 
symmetrical functions are represented by graphs where the 
number of nodes grows at most at a rate proportional to the 
Square of the number of arguments (n’*). The BDD therefore has 
a more compact representation. Secondly, the reduced form of 
BDD 1S canonical. (i.e., every function has a unique 
representation), [2]. 

The binary tree can be explained in the form cf a Shannon 


expansion: 
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Fo 
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Figure 1. Shannon expansion. 











Using the logic circuit of Figure 1 to implement the 


expansion, we have 


In the same way, we can expand the new subfunctions as follow: 


FOO SS EangemOR cha fz 
FOl = X3'f919-XOR.X3f5,, 
£10 = X3'f459-XOR.X3f,9, 
fll = x3,'1f,,)-XOR.x3f,,, 
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Figure 2. Representation of logic function using Shannon 


expansion. 











If we replace each node by a 2 input multiplexer, we have the 


following network: 


X1 











fOO0OO §=©foo10 «701900 §=—f0110 #1000 = f010—S ss f1100 ~~ #1110 
{0001 f0011 f0101 f0111 1001 f1011 1401 1111 


Figure 3. Realization of logic function using multiplexers. 


A BDD is a way to represent a given function using a 
binary tree. A binary tree has a root node which is unique in 
that it is not a child node of any other node. The root node 
of a BDD represents the given function f£(x,,x5,...x,). The 
left child of the root represents the  subfunction 
£(0,X5,X3..-.*,) and the right child f(1,x,,x,...x,). Similarly, 
the grandchildren, great-grandchildren etc. represent the 
subfunctions associated with x, and x, and so on, until all 
variables are assigned. The leaf nodes, which have no children 
are assigned with constant O and 1. 


ee 


AND OR 
Figure 4. The BDD for a AND and OR’ function. 


Figure 4 shows examples of the BDD for the AND and the 
OR function. In this case, the ordering of the variables 
makes no difference in the BDD structure. Ref [2] showed an 
interesting example of a BDD which is highly dependent on 
ordering of the variables, (See Figure 5) The functions x,x, 
+ X3X, + XoXg and xX,X, + XoXo + X3X_ differ from each other 
only by a permutation of variables, yet the first BDD has 8 
nodes while the second has 16 nodes. From this example, one 
can see that a poor choice of ordering can have very 


undesirable results. 











0 


Figure 5. Example of argument ordering dependence. 


Besides the difference in complexity in terms on number 
of nodes, the functions x,x, + xXx, + X;xX, has crossings, while 
the other does not. This shows that BDDs without crossing can 
be found by varying the order of the variable, if they exists. 

In the LSI implementation of a network, crossings are 
expensive because they require additional channels and 
increase delay. In submicron LSI, the delay in the 
interconnections are comparable to the delay in the logic 
elements. For a complex network, the delay caused by 
crossings can be unacceptable. It is therefore understandable 
that one of the important goals in circuit design involving 
BDD’s is to eliminate, or at least reduce crossings in the 
network. 

Sasao & Butler [4] define a restricted planar BDD: 





Definition 1: A function has a restricted planar BDD if 
there exists a BDD without crossings, where the 1 edge emerges 
to the right of the node and 0 edge emerges to the left of the 
node; and the constant 1 is in the right of the constant 0. 


This definition will be used to describe a planar BDD in 
this text. Until recently, only BDDs up to 5 variables have 
been shown to have no crossing [4]. The progress in this 
field is slow because manual enumerations of BDD are very 
tedious. To facilitate further research related to BDDs, the 
process of drawing BDD, searching for those without crossings, 
and producing their orderings, counting of their nodes etc. 
must be computerized. 

The development of such a tool is one of the main aims of 
this thesis. The application of this tool is demonstrated by 
characterizing a class of functions. The following 
definitions in accordance with [1] are used in the subsequent 


discussion: 


Definition 2: A switching function f£f(x1,x2, ....xn) is 
a mapping of f£:B"->B, where B = {0,1}. 


Definition 3: A threshold function f£(x1,x2, ....xn), has 
the property that f£ = 1 iff WX, tWy-1X,-y+----+W,x, > T, where 
T and w, are integers and the logic values, 0 and 1 of x, are 


viewed as integers. 


The value of WX) t+Wy-1X,-1+----+W ,xX,, for some assignment 
of values x1,x2, ....and xn is called the weighted sum. A 
threshold function is completely specified by a weight- 
threshold vector, (w,,W,.4,--..,W, ;T). AND and OR are special 
examples of threshold functions. Since any Switching function 
can be realized by a combination of AND, OR and inverter, any 
Switching function can be realized by a network of threshold 











functions. In threshold logic, in general, each logic gate 
represents more than AND,OR or NOT does in conventional 
Switching theory. Accordingly, the number of gate that 
realize a given function is often less than the number of 
gates required by the AND,OR, NOT or other conventional gates. 








II. REQUIREMENTS OF THE ANALYSIS TOOLS 
This thesis describes a tool to analyze BDD’s for 
threshold functions. It can be extended to evaluate other 
functions. The tool is designed with the following 


considerations in mind: 


1. Inputs 
The parameters describing the threshold functions such as 
the vector (wW,,Wy.1 ----,W2,W,; T) will be requested once the 


program is started. Several options are provided to tailor 
the analysis to certain requirements. Since some analysis 
requires the threshold to be varied, the option to vary the 
threshold from some minimal value to a maximum value is also 
provided. If the aim is to search for planar BDDs, the 
option for producing various ordering of the variables in 
lexicographical order is available. The option for drawing 
the BDD and printing in a text file can be activated to help 
in visual analysis of BDD’s. However, for certain BDD 
structures, when searching a threshold function of many 
variables, the number of planar BDDs can be very large. A 
preview of which BDD is planar can be made by requesting a 
listing of planar BDDs first. 

2. Processing 

The program is designed to evaluate a threshold 
function and provide the result in the form of a truth table 
for further BDD analysis. For other types of functions, the 
routine to convert the particular function to a truth table 
must be written. The truth table can be a large table 
having 2” entries, where n is the number of variables. Truth 
tables are useful for checking or interfacing with other 
programs. It is also used in the generation of the initial 
BDD. 

The program simplifies the initial BDD to a form called 
reduced ordered BDD. To allow drawing of the BDD, the reduced 
BDD is laid out in the form of a binary tree. 











At present, the program is able to analyze BDDs of up to 
9 variables. If threshold functions of more variables are 
needed, the program can be expanded. The program is designed 
to detect any crossings in the BDD. This feature would help 
in searching for planar BDDs in any type of function, and is 
useful for the design of fast logic circuit. An important 
output is the number of nodes of a reduced order BDD. This 
feature is useful since the compactness of a BDD is given by 
the number of nodes in the BDD. 

The program produces all the possible orderings of 
variables to search for planar BDD, if it exists at all. As 
discussed before, some orderings yield planar BDDs. 

3. Output 

The program provides the following output for BDD 


analysis: 


® Layout for drawing of the BDD 
@® Number of nodes for each BDD 


® List of BDD’s without crossings, showing the ordering 
for each BDD 
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IIIT. ALGORITHM & IMPLEMENTATION IN C++ 

The algorithm to enumerate binary decision diagrams and 
analyzing them is described below: 

1. Generate Lexicographical Ordering 

This module is required for the enumeration of BDDs. It 
produces all orderings of the variables. If two variables 
have the same weight, only the first ordering is used; all 
subsequent orderings of variables corresponding to this same 
arrangement of weights are ignored. For example, with weight- 
threshold vector (3,2,1,1; 5), ordering x,x,x.x, is ignored, 
and only ordering x,x3x,x, is used. The program run time can 
be very long as the number of variables increases since there 
are n! threshold functions to enumerate. 

2. Enumeration of Threshold Function 

An important part of the program is to generate a 
threshold function with weight-threshold vector 
( Wy,---,Wo,W,;T ) which is stored in an array to be used for 
binary tree generation at a later stage. 

3. Generation of a Binary Tree 

A binary tree is generated with 2'*1) nodes. Each of the 
nodes of the binary tree has the following properties: 


® Index: Each node is assigned an index that increases 
from top to bottom and from left to right. This 
identification is used later to check whether a node 
can be merged. Figure 6 shows an example of the 
indexing of the nodes for a binary tree that describes 
a 4 variables function. For every node of index i, its 
left branch has an index of 2i, and its right branch 
has an index of 2i+1. 


® Pointer to left branch: Pointer that point to the left 
node immediately following that node concerned. 


@® Pointer to right branch: Pointer that point to the 
right node immediately following that node concerned. 


® Data: A data value, DATA, associated with all nodes 


other than leaf nodes (those with index from 2" to 
2 (n+1) ) iS initialized with the character "*". In the 
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Subsequent steps, other values such as 0,1 or an 
alphabet "a" ..."j" will be assigned to indicate the 
Size of the subfunction that has been merged. A leaf 
node corresponds to some assignment of values to all 
variables. Each such node is labeled by the value of 
the function for that assignment. | 













5 f01 f10 fi 


incre 


fO000 0001 f0010 0011 f0100 10101 10110 0111 f1000 #1001 f1010 (1011 1100 f1101 frt10 F144; 





Figure 6. Indexing of Binary tree. 
4. Simplification of BDD 
Figure 6 is also an example of a full binary decision 
diagram of a function. Such a BDD has many redundant nodes 
and can be simplified. The algorithm to do this is as follow: 


a. If the value of a subfunction, which is 
stored in DATA, is a constant 0 or 1, terminate 
that branch. The DATA of the nodes of the binary 
tree that proceed after the constant 1 or 0 will be 
filled with blanks to show that the corresponding 
branches are terminated. 
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b. If the current subfunction f, is the same as 
one already generated f. , move the current branch 
f; (£;) over to the other branch f,. For example, if 
fOl is the same as £10, move £10 branch over thus 
merging the 2 nodes. The BDD resulting from the 
above simplification is called a Reduced Order 
Binary Decision Diagram (ROBDD) [3]. Letters "a" 

"7" are used to indicate the size of the 
subfunction merged. If the merging is done for a 
subfunction with only two children and 0no 
grandchildren, then the merged subfunction will be 
assigned a character "a". Similarly, for nodes with 
two children and 4 grandchildren will be assigned 
the character "b" if it is a candidate for merging. 
etc. Such a scheme is implemented for nodes that 
have up to 2? offsprings. See the example in table 
1. The merging of the subfunctions is done from 
"bottom up", which means that comparison of 
subfunctions are done first for subfunctions with 3 
nodes. All subfunctions with 3 nodes are compared 
and merged if found to be identical. A merged 
subfunctions with three nodes will be assigned with 
Character "a". After that the next level of 
subfunctions with 7 nodes will be compared. The 
procesS goes on until all subfunctions’ are 
evaluated for their potential to merge. The size 
of the largest subfunction compared is 241. 


Number of nodes in 
subfunction merged 





ac 


Table 1: Example of descriptor for subfunction. 


3 





oP If two subfunctions of one node n are the 
Same, i.e. the node’s left branch ny and the right 

branch n, are the same, then extend the branch to 
n down to the ny and n, eliminating n. (i.e., do not 
use a multiplexer ). 


5. Counting Nodes 

The number of nodes is a useful measure of the 
compactness of a BDD. C++ has a recursive binary tree 
traversal function that is very convenient for this purpose. 


A valid node has the characteristics of DATA = "*" with two 


children branches. 
6. Check Crossings 
inorder traversal is used in the checking of crossing. 


The steps for inorder traversal are : 


Traverse left subtree using an inorder traversal. 


Visit the root node. 
Traverse right subtree using an inorder traversal. 
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Figure 7. Example of inorder traversal. The order of 
traversal for the above Binary tree is 8-4-9-2-10-5-11-1-3-14- 
TSi5% 


Using the recursive C++ function "inorder traversal", 
the DATA elements of each node of the Reduced Ordered Binary 
Diagram are checked consecutively for a transition of "0" 
to "1" . Not all nodes have a 0 or 1. When aoO or 1 is 
encountered, it is recorded to check for 0-1 transition. For 
a reduced order planar BDD, there should only be one such 
transition during the complete traversal of the binary tree. 
i.e., All nodes preceding the last "0" must be "0", "*" or any 
symbol representing a merged subfunction. Likewise, any node 
after the first occurrence of a "1" must be either "1" or "*" 
or any of the symbol representing a merged subfunction. A BDD 
with more than one 0-1 transition in a inorder traversal of 
the reduced binary tree will exhibit crossings, since all 
leaves with "0" can only be joined to each other with 
crossing, and likewise those leaves with "1". Figures 8 
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Figure 8. (Step 1.) Enumeration of the threshold function 


with weight-threshold vector (5,3,2,1,1;7). 
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Figure 9. ( Step 2.) Simplification of BDD. 





0 1 


Figure 10. (Step 3.) Printing of simplified BDD layout. 
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X1 


X2 = & 


X3 CY & 
X4 Ee 


X5 
0 1 
Figure 11. (Step 4.) Completion of Planar BDD 
The current implementation of the algorithm is able to 


analyze threshold functions of up to 9 variables. Figures 12- 
15. are examples of 6,7,8 and 9 variable Fibonacci threshold 


functions represented in the form of planar BDD. 
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Figure 12. Example of a planar BDD for function with weight- 
Ghreshold vector of (8,.5,.3,2,/i,i13 12). 
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Figure 13. Example of a planar BDD for threshold function with 
weight-threshold vector of (13,8,5,3,2,1,1; 23). 
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Figure 14. Planar BDD for function with weight-threshold 
VECEOrD (21135 6):5)15;2744 > 3S). 


eae 





Figure 15. Planar BDD for Fibonacci threshold function with 
weloht-threshold vector (34,21,13,875,3,2;i,1; 350). 
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Figure 16. Planar BDD for threshold function with weight- 
threshold vector .6,7557-4,4;,3,3,2 2: 8). 
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7. Output 

There are two main output forms: 

a. Binary Decision Diagram: The nodes of the planar BDD 
is laid out in the form of a binary tree in ASCII 
text. For a large tree (i.e. of 7 or more variable 
threshold function), the tree should be viewed or 
printed out in compressed font due to the limited space. 
The nodes can be joined together in the same way as a 
binary tree. The merged subfunctions are represented by 
alphabets a,b,....j. ‘’a’ corresponds to the subfunction 
with two leaf nodes, b correspond to a subfunction with 
two children nodes, each with two leaf nodes etc. These 
nodes marked by a letter (a,b,...,j) should be joined to 
the equivalent subfunction that it represent on the left 
of the binary tree. All nodes with leaf nodes of ‘0’ and 
‘1’ should be joined to sink nodes ‘0’ and ‘1’ at the 


lowest layer of the binary tree. 


b. Listing: The listing of the planar BDDs without 
crossing are printed to a ASCII text file in the same 
directory as the program. Other characterization, such 
as the number of nodes, ordering etc. are also listed. 
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Iv. ANALYSIS 
A. PLANAR BDD FOR CLASSICAL THRESHOLD FUNCTIONS 


Muroga [6] defines canonical function as a function of n 
variables with w, > Wy >...>w,. He lists all canonical 
positive threshold functions with five or fewer variables. 

For six variables, only canonical positive self dual 
functions are shown. A dual of the function f(x) is defined 
as 

£%(x) = £' (x’) 
where f’ is the complement of f and x’ = (x,', x,’,...x,’) and 
A function is self dual if 

£o(x) = £ (x) 
x,’ is the complement of x,. 

Using the methodology described in this thesis, the 
ordering that produces planar BDDs for each of the listed 
threshold functions is found. The orderings that correspond 
to these planar BDD are listed in Appendix (A). From this, we 
can state 


Theorem 1 : All threshold functions of up to five 
variables and all six variables canonical positive self dual 
threshold functions have a planar BDD. 


Based on the ordering suggested by the program written, 
the BDD can be either derived manually or by means of the BDD 
program whose source code is listed in APPENDIX (C), for all 
threshold functions listed in the APPENDIX (A). This table 
greatly simplifies the design of efficient logic network 
based on planar BDD for common threshold functions. Giving 
the number of nodes that corresponds to each planar BDD helps 
in further optimizing the logic design, by indicating which 
ordering has. the minimal number of nodes. 
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Theorem 2. All threshold functions of up to five 
variables and all six variables canonical positive self dual 
threshold function have a planar BDD with the minimum number 


of nodes. 


This shows that planar BDDs are amongst the most 
efficient representations of the threshold function with a 
small number of variables. This is shown to be the case for 
threshold functions up to 8 variables. 

The percentage of planar BDD decreases rapidly as the 
number of arguments increases, which means that planar BDD are 
extremely rare as the number of arguments increase. The 
orderings for planar BDD cannot be easily derived 
heuristically for threshold functions of many variables. 
There is often more than one minimal BDD. The following table 
Shows the frequency of planar BDD versus all possible BDD’s 


Average no. of BDD’s 
that are planar 


Average no. of BDD’s 
that are planar and 


minimal 


Percentage planar 
BDD ({ 3% ) 





Table 2 : Planar BDD is very rare as the number of variables 
increases. Column 6 refers to 6 variable canonical positive 
self dual threshold functions. 
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The total number of permutations refers to the total 
number of permutations on on variables that produce unique 
functions, for all n-variables threshold functions listed in 
Appendix A. The average number of permutations that are 
planar is found by dividing the total number of planar BDDs 
by the number of threshold functions. When such BDD has 
minimum number of nodes, it is classified as planar and 


minimal. 


Percentage of BDD that are Planar 
100 


Percentage of Planar 8DDs 
gi 
© 


1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 
Number of variables 


Figure 17. Percentage of planar BDD in all unique permutations. 


zy 








Figure 17 is a graph of percentage of permutations that yields 
planar BDD versus the number of variables. Planar BDDs for 
threshold functions of more than 8 variables are very rare. 
Although exhaustive search planar BDDs for all threshold 
function with 8 variables is not done, search on a_ small 
sample of 8 variable BDD, such as the one shown in Figure 16 


shows that percentage of planar BDDs approaches zero. 
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B. FIBONACCI FUNCTIONS 


Definition 4: A Fibonacci number F; is specified by a 


recursive relationship: F, = Fi, + Fi., , with F,=Fo=1. 


Definition 5: A Fibonacci function is a threshold 
function with weight-threshold vector (F,,F,.1,....F2,F,:T), 
where F, is the i-th Fibonacci number and 0 < T < Feo? 


A Fibonacci function is a threshold function with weight- 
threshold vector Fi ea eeah oy hye) y where F; is the i-th 
Fibonacci number and 0 < T < F,,.. The following have been 
observed for Fibonacci threshold functions of up to 9 
variables. 


Theorem 3. All Fibonacci threshold functions with 
thresholds of one to the largest threshold and with an 
ascending ordering of variable have planar BDD.s 


This is also stated in [1]. The largest threshold is 
the same as the largest weighted sum, 1+1+2+3+....+F, =F, 5. 
Table 1 in Appendix B shows the number of nodes of Fibonacci 


functions for threshold of 1 to the largest threshold. 


Theorem 4. Amongst the orderings that yield the smallest 
number of nodes are (1.) descending and (2.) ascending order 
of weight. 


When the threshold is 1, the ordering does not matter at 


all, since this function corresponds to a OR function. When 
the threshold is at its maximum value, the planar BDD 
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correspond to a flipped image of the BDD with threshold of 1. 
It can be seen from the listing in Appendix B that these 


always correspond to the minimal BDD. 


Theorem 5. The compactness profile of the planar 
Fibonacci BDD’s with thresholds from 1 to its largest 


threshold is symmetrical. 


This symmetry can be seen clearly graphically as shown in 
Appendix B. Figure 18 shows an example of Fibonacci functions 
of 9 variables, with the plot of the number of nodes versus 
the threshold value of 1 to their maximum. The distance from 
one minimum to the next one in the compactness profile 


exhibits an interesting characteristic. i.e., 


Distance of it? trough to (i+1)® trough 
~ Golden ratio 





Distance of (i421) th trough to (i+2) th trough 


Fu 
with Golden ratio = lim ————— = 1.618 . 
n->@® Fn-1 


where F, is a Fibonacci number . 
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Figure 18. 
weight-threshold vector of 
thresholds (T) of 1 to 54. 
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Figure 19. Distribution of Fibonacci function by Nodes and 


Variables. 


Figure 19 shows the distribution of nodes in the BDD’s of 
Fibonacci functions as enumerated by the BDD program. The 
number of variables and the number of nodes in the BDD’s are 
plotted horizontally, while the number of Fibonacci functions 


is plotted vertically. 
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V. CONCLUSION 

An effective method of finding planar BDDs for threshold 
functions is developed. The algorithm is implemented for 
threshold functions of having up to 9 variables. The source 
code in Borland C++ is shown in Appendix C. It ais 
demonstrated that Fibonacci threshold functions having up to 
9 variables have planar BDDs. The Fibonacci threshold 
function is also characterized using the algorithm. It is 
also found that all threshold functions of having up to 6 
variables listed in [6] have planar BDDs. The ordering that 
produce planar BDDs are catalogued in Appendix A for easy 
reference. 

Follow up research in this area may produce an algorithm 
to convert a BDD to its sum-of-product expression. This would 
facilitate research in more effective representation of 
Switching functions. Currently, the program does not accept 
other functions as input. It can be made more versatile if 
input in the form of a sum-of-product expression to the above 


described program can be implemented. 
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APPENDIX A- ORDERING FOR PLANAR BINARY DECISION DIAGRAM 


The following tables, extracted from [6], list all 
threshold functions of 3 to 6 variables. It is found that all 
the threshold function have planar BDD. Many of the orderings 
are symmetrical due to identical weights. Only one of the 
symmetrical BDD is listed here. In the following table, in 
the ordering of 321 means x,x,x,, with x, as the root node, and 
the xX, aS sink node in the BDD. For each threshold function in 
the table, a disjunctive form is expressed with the 
subscripts of the variables; for example, 1v23 denotes 
x,.OR. (x,.AND.x,). 


The planar BDD for the following three variable threshold 


functions are found and tabulated with their orderings and the 
number of nodes, which indicates their compactness. 


S/N | Weights; Ordering | Number of Threshold Functions 
Threshold nodes 
Ea Ce eee 


BW a Be 











S/N | Weights; Ordering | Number of Threshold Functions 
Threshold nodes 
dle 221125 4321 123Vi24 
2431 
1234 
2: 224172 4327 lv2v34 
2435 
2143 
ca LiL. 3 any 123v124v134v234 
order 


gs a I es any Ll2v1l3v23vl14v34v24 
order 
oe ZL1 134 4321 7 EZ Svigayi 34 
1432 7 
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2134 
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en 


L234: | otis 3 4321 l2vi3v1l4v234 
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12: | ek 4321 1l2v1l3v14 
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1324 


Table 2. Planar BDD for four variables threshold functions 
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Table 3. Five variables threshold function. 
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Table 4. Planar BDD for 6 variables threshold functions 
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APPENDIX B. COMPACTNESS OF BDD FOR FIBONACCI FUNCTION 


The following table illustrates the compactness of the 
BDD for fibonacci threshold function of 9,8,7,6,5 variables, 
which corresponds to column B,C,D,E,F in the table below. The 
profile of the number of nodes is symmetrical as the threshold 
varies from 1 to its maximum. The plot of theses nodes shows 
that the ratio of one gap to the next gap approximates the 
golden ratio. (Gap = distance from one minimal point to the 
next ) 


7 
| 
a 
9 
7 
9 
9 
5 
9 
qd 
7 
7 





Table la. Number of nodes of fibonacci function at 
different thresholds. ( To be continued in table 1b. ) 
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Table 1b 
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Conta 
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Table lc. ( Cont’d ) Number of nodes of fibonacci 
function at different thresholds. 


Figure 1. to Figure 5. show the plots of the profiles for 
the number of nodes versus the different thresholds of the 
Fibonacci Tunctron. 
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Figure 1. Compoactness protiie for F(5;3,2,2;12T) 
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Figure 2. Compactness profile for F(8,5,3,2,1,1;T) 
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Figure 3. Compactness profile for F(13,8,5,3,2,1,1:;:T) 
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8 variable Fibonacci Threshold Function 
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Figure 4. Compactness profile for F(21,13,8,5,3,2,1,1;T) 
9 variable Fibonacci Threshoid Function 
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Figure 5. Compactness profile for F(34,21,13,8,5,3,2,1,1;T) 
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APPENDIX C. SOURCE CODE FOR BDD PROGRAM 


[EERE KRKKKEERERERKE KEKE EEE RR EE EEE EEK EE HEE KEKE KKK KKK KR E EKER / 


is Filename: TREE.H ‘di 
j> é 


[KKK KERR KKK KK KK RK RRR RRR RRR RRR RR RR KK RR KR KKK KK KKK KK / 


#define NULL 0 
typedef char DATA; 


struct node 


{ 


DATA ad; 

int index; 

Struct node *Lere; 
struct node *Prgnts 


tie 


typedef struct node NODE; 
typedef NODE *BTREE; 


[RREKKKKKKKKKKKKKKKKEK Prototypes FA TE A I RIE NG So eM Te ee ee 


void print_to_file(int node_remain, DATA p[], BTREE root, int 
numvar, char filename[], int threshold, int wf], int 
x_order[]); 


BTREE init_node(DATA dl, int index, BTREE pl, BTREE p2); 
BTREE new_node(voidq) ; 

BTREE create_tree(DATA a[], int i, int size); 

void inorder(BTREE root, int *cross_flag); 

void preorder(BTREE root, int *node_remain); 

void postorder (DATA p[],BTREE root,int *node_cnt, int num_var 
e 

DATA merge (BTREE node_ptr, int *node_remain) ; 

void print_tree(DATA p[], BTREE root, int num_var, int 
x_order[]); 

void cross_test (DATA d, int *cross_flag_ptr); 

void merge_check( DATA *f_1[], int num_var); 

DATA node(DATA node_loc[], BTREE root): 

void data_set(int tf, int *threshold, int *woinl])- 

int sfl_cmp(DATA sfl1[]); 


9. 











PERE ALLER EES EARS RR RE RRR RT ROR AR AR RD RR SRE A I Ie eR a IK 


/* Filename: BDDANATL.C as 
i a7 
/* BDD Tree Analysis Program for < 10 Variables ae: 
is ey 
/* Pounctions: my 
/* a.Permutate all orderings of a Threshold Functions “es 
/* b. Evaluate the Threshold Functions a 
/* c. Generate BDD for each Function. oe 
/*= Os Simplify the BpD: uy A 
/* e. Examine for Crossings. “7 
/* d. Store results in ASCII file. ia 
as The results consist of * / 
eal wy é 
dial ae Ordering of variables that produce planar BDD x / 
<> 11. Compactness of BDD in terms of no. of nodes ug 
ae 111. Binary Decision Diagram of the Function ry 


LEER TREEAAREKE ERE EAR ELLE RALAA ERR CELE RRRRAR E AERA R A AR A RY 


#include 
#include 
#include 
#include 
#include 


void main 


{ 


Int 


iy 


cha 


int 


r 


©, 


"tree.h" 
<stdio.h> 
<stdlib..A> 
<math.h> 
<string.h> 


min_node=1000, 


Carry [20 |, 


cmp_vector[20], 


found_flag, 


i. 
a 


1, 


f 


Mm, 


5 


sf1[2050], 


Si2(2050] ; 


f_ size; 


/* Counter for searching for minimum 
no. of nodes */ 
/* Used in evaluation of threshold 


function ae 4 

/* Vector used for comparison of 
sub-functions ae 4 

/* Flag. to andicate. that .a “valid 
ordering is found*/ 

/* Index */ 

/* Index */ 

/* Index */ 


/* Size of subfunction to be compared 


J String describing FLYSc 
subfunction to be compared */ 
eS String describing second 


subfunction to be compared */ 


/* Index */ 


56 








DATA 


FILE 


BTREE 





w[20], 
node_remain=0, 


*node_cnt, 
num_var, 
cross_flag=0, 
xorder|( 201, 
threshold, 
bin_wt [10], 


bin_index[2050], 
Cmpstr{2050 J: 


x10] [L024 ) » 
f_vector[2050], 
IMpUt Ok 
filename[8], 
Sl 20SC) ; 
PLO; Z0S0), 


p[2050]; 


“OLD; 


root; 





/* Index */ 
/* Index */ 
f*™. THGGx -F/ 


/* Enumerated threshold function as 
weighted sum */ 

/* Weights of threshold function 
input */ 

/* Weights of threshold functions */ 
/* Counter for number of nodes in 
binary tree */ 

/* Counter for number of nodes left 
in binary tree */ 

/* Number of variable in a threshold 
function */ 

/* Indicate that there is a crossing 
if set */ 

/* Ordering of the weights in a 
threshold function */ 

/* Threshold in a threshold function 
vy 

/* Binary weights used for computing 
leaf node address */ 

/* Leaf node address * / 

/* Temporary string used in BDD 
manipulation */ 

/* Max number of variable =10 */ 


/* Vector describing threshold 
function to be converted into binary 
tree*/ 

/™ Input. control tiag a 

/* Output filename */ 


/* Threshold function value in 1s & 
Os,merging are marked here */ 

/* Unchanged Original Threshold 
function vector */ 

/* Data describing a node in binary 
tree */ 


/* Pointer Lor output file */ 


/* Root address of binary tree */ 


oi | 


PE pee Data 20pUL, Sees Sen ees wes 


print£(" \n Number of variables : "); 
scanf£(" $da",&num_var) ; 


for(1=0; i<num_var; i++) 

{ 

printf("Weights :w%d=",1); 
scanf("%d", &w_inf{i]) 


} 


printft(" \n Enter threshold ="); 
scanf("%d",&threshold) ; 


for(1=0; 1<num_var; i++) 


{ 

printf("Order X:w%td= ",1); 
scanf("%d", &x_order[i]); 

} 

printf("Threshold Function : F( "); 
printf (" Weight : "); 

for ( 1=0; i<num_var; i++) 
Prinli(" sa: “ won) 

prantt(" +7. $d ).",threshold): 

Princt ("ns 

brinti(* Order 22) 

for ( 1=0; i<num_var; 
printf("%d ",x_order[i]) 

Preintt( \0"); 

/* Loop for various threshold ,provide the 


thresholds that should be evaluated with the 
functions here */ 


for (threshold=1; threshold<13; threshold++) 
{ 


ofp=fopen("#node.5ve", "a"); 


fprintf(ofp, "Fibonacci function with %d variables 


threshold = $d ) \n",num_var,threshold ); 
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range of 
febonacci 


, and 





for (1=0; 1<20; i++) 
: 
carry [1]=0; 
1f (i<num_var) 
x_order[i]=num_var-i; 
else 
x _order[i]=0; 


x_order[0]=num_var-1; 
/**** generate increasing number for permutation generation */ 


do { 
x_order[0]++; 
for (1=0; i<num_var; i++) 
{ 
x_order [i]=x_order[i]+carry[1i1]; 
Carry [1]=0; 
if (x_order[i]>num_var) 


i 
i ie 


x order [i] 
carry [i+1] 
} 


} 
PPR SCarL se “*** Tier oul. he. permitcatrton: ert er/ 
for (j=l; j<=num_var; j++) /* number to compare */ 


{ 
found_flag=0; 
for (1=0; i<num_var; i++) 
{ 
if (x_order[i]==)j) 
{ 
found_flag=1; 
break; 


} 


if ‘(£ound flags=0) es Lf found_flag=1 after 
comparing 
all elements of the array, it 
1s a valid permutation */ 
break; 


ope, 














1f (found_flag==1) 
{ 


/* Here a permutation of an ordering is computated. The next 
thing to be done is to find out which of the Fibonacci 
function in this ordering are planar, and has the minimum 
number of nodes at the same time */ 


min_node=1000; 


Draner(*\n\n*) > 
printf("\n Threshold = $d Ordering = ",threshold 3 
for (j7=0+ Jenumo vars 344) 
print£("%sd",x_order[num_var-j-1]); 
PElner (™\a\n™): 


[xxxxekekEKKkKDermutation ends here PE I Ae CR IER I eRe OI He 


FOr (2=0; 22num vars 244) 
w([1)swlinl| (x.order[iJ=-1)] > 


/*---— Initialisation ~-+------------+-----+---------+.4 iy 4 


Lor (220s. 2210s 243) 
bin_wt [num_var-i-1]=pow(2,i); 


for (2203 2= (pow (2,20) ) + i44) 

bin_index[i]<=0; 

tor (ie0y aenum vars 344) 

{ for (j=0; jJ<pow(2,num_var); j++) 
Say aCe 


For (120 A<num-vare 244) 

{ for (j=0; j<pow(2,num_var); j++) 
Sa eos 

} 


for (j=0; j<pow(2,num_var+1); j++) 


EL )e05 

Pigl=" 2% /* initialise all nodes with the symbol 
blanks */ 

f1[3]=EOF; 


£f_lo[j3]=EOF; 
} 
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Dizo4cie"\0": 


f Cae eae Generate X0 Xl... Kn where n=num_var ---------- | 


FOr (1=0)) 2enum vars 4-64) 


EOr(K=0; ke (pow(2,2) )s kee 4 
{ 

> he 
J=J+1; 


} 


ror (k=0; ke<(pow(2,3))* kee 3 
{ 

> Nd as 

J=j+1; 

} 


} while ( j<pow(2,num_var) ); /* while loop */ 
y Ee £Ot Loom: ey 


art aan generate results of threshold function ------- mh 
for ( 3=0; j<pow(2,num_var); j++) 

FOr (i=0f a= num vars 144) 

{ 


£(3]=w[i]*x[ (num_var-x_order[i]})][j]+f£[3]; 


bin_index[j]=bin_wt [i] *x[ (num_var-x_order[i])] [j]+bin_index|[ 
113 


} 
for { 3=07 4< (powt2 num var) )3 j++) 
{ 
if ((f[j] == threshold) | | ( f£{[j] > threshold@) ) 
{ 
f_1[(bin_index[j])]='1':; 
f_lof{ (bin_index[j])]='1’; 
rea is ie 
} 
else 
{ 
f_l{(bin_index[j])]='0'; 
f_lo[ (bin_index[j])]='0'; 
st Basle ee 
} 
} 
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/* Vector of threshold functions value is formed, the next 
step 1s to simplify by comparing adjacent subfunctions. 
Identical subfunction will be combined. */ 


fOr (j=l ‘enum vars 744) 


{ 


for (1=0; 2< (pow (2 -num var) ) = 144) 


{ 
Sra] S407 
StZ [aa " V0": 
} 


sf_size=pow(2,j); 
1=0; 
while (1<(pow(2,num_var) )) 
{ 


for (k=1; k<(sf_size+i); k++) 


{ 
sf1l[(k-i] = f_lo[k]; 
sf2[(k-1] = f_lo[(k+sf_size)]; 
} 
if ((stremp(sfl,sf2) == )&& (stremp(sf1,"00") !=0)&&(stremp(sfl, 


"Q0000") !=0)&&(stremp(sf1,"00000000") !=0) &&(stremp(sf1,"00000 
00000000000") !=0)&&(stremp(sf1,"0000000000000000000000000000 
0000") !=0)&&(stremp(sf1, "00000000000000000000000000000000000 
00000000000000000000000000000") !=0) && (stremp(sf1,"11") !=0) && 
(stremp(sfi,"1111") !=0)&&(stremp(sf1,"11111111") !=0)&&(strem 
Distl (Titi Lia 0 ee (Sbremm (stl, "Titi tii rid 
eee PT Lay be0V ee Stremp (Sb i 114 9 ed 
pt Ee ee 0g 4 


{ 
1f ((sf1!="00000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000") && (s£1!="1111111111111111111111111111111 
ee Ee eel Oa ct De Bp tS tO Zs Ue Yas i Ps DD Ea a hp Ks Ui a se Pet Oe tc 
ss pW dG oD Ut ZU et ps Ut Wt et et Bp Ui sa 2 Dae TU UU Ds TR 
{ 

1£((sfl!="00000000000000000000000000000000000000000000000000 
DODDDD00000000000D000ND00000DND0ND0ND00000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
00000000000000000000000000") && (sf1!="11111111111111111111111 
i 3 0 es Ut Ot a Ee Bt FU Ug Fo 0 et pW Eg Bs Ds UD Ut Ft es UT Us UD Ws Es Ts oD 
is BW 5 et Fst Ds i Wt es Te et DoD Us Pe i Rit eS ps Ups Ee et eV UW Dis Us Ue es CU 
i Me Be ele Di Eo 2 De Cs 6 BK is Ee Uti a Et ig pF Ts eV Os Ws ya Ee Ui Ys Ep Ts Wt sO GE es 
US Bt et Ut CB Oe Det Ke UE as et se pW Di Ut Ve Ts eB et TD ig EU Ue Ut Bs oe 


{ 
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for (m=(1+sf_size); m<(1i+2*sf_size); m++) 
{ 

1f (sf_size==2) 
Eo | i Sen 

else if (sf_size==4) 
sme OB 0 ea oe 

else if (sf_size==8) 
j eee OB oc eae oa 

else if (sf_size==16) 
pe ea ea Ba 

else 1f (sf_size==32) 
t.1.im)="e** 

else if (sf_size==64) 
Elm) eS 

else if (sf_size==128) 
time" os 

else if (sf_size==256) 
f_l{m]=‘h’; 

else if (sf_size==512) 
f lf{mj='1’; 

else 1f (sf_size==1024) 
same 0 Fee cae 

} 


} 
} 
} 
1=1+2*sf size; 
} 
} 


f_l[pow(2,num_var)]='\0’'; 


fOr: (jels jenumvare 744) 


{ 


for (1=0; 1<(pow(2,num_var)); i++) 

{ Sti lij="\0": /* Initialisation 
sf2[iJj=’\0'; 

} 


sf_size=pow(2,j); 


1=sf_size; 
while (i<(pow(2,num_var) ) ) 
{ 


for (k=1; k<(st£..s1Z641) > kes) 


{ 
sf1[k-1] 
sf2[k-1i] 


f lo[k]; 
f_lo[ (k+sf_size)]; 


} 
ye de Sei emp (st 1): 7 
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if 








if ((stremp(sf1,sf2)==0)&&(stremp(sf1,"00") !=0)&&(stremp (sf, 
"0000") !=0)&&(stremp(sf1,"00000000") !=0)&& (strcmp (sf1, "00000 
00000000000") !=0)&&(stremp(sfl1, "0000000000000000000000000000 
0000") !=0)&& (stremp(sf1, "00000000000000000000000000000000000 
00000000000000000000000000000") !=0) && (strcmp (sf1,"11") !=0) && 
(stremp(sf1,"1111") !=0)&&(stremp(sf1,"11111111") !=0) &&(strem 
D(stiy LITT Ate i )t0) ee (st remote to 11 
PTTL 0) ee (Stremo (Stl. 81 14 aa Pe 
i Ht es Ut Le Dt es Dp Us Tt Kt et Ka et sD DDO eo WRC en 


{ 


1f((s£1!="00000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000") && (sf1!="1111111111111111111111111111111 
Es Et i tt eB Di ep Ut Et Dp he Wigs Di 9s Kit tO WD pV oe het eT tO 
BE sb it pt te Et aE Dt at eK Ep CB nD Ut Che Us Et Ti Ft tt WF Ts | 

{ 


1f((sf1!="00000000000000000000000000000000000000000000 
DO0DD00000000000000D0D0DDNDDNDND0DNDDDNDD0000000000000000000000000000 
JDODDD0D0DDND0D0D00000NDDNNDDNNDDDNDD0DD0000000000000000000000000000 
0O0DD00000000000D0DNNDNDNDNDDNDNDDDDN0ND0000000000000000000000000000 
00000000000000000000000000000000") && (sf1!="11111111111111111 
LL4141421111121111111111111121111111111111111111111111111111 
ites a at tl sO es a Ds Ut Fp Ts Es i Ts ig ot Ts ps ff Ke i TV it tet EDC it PD cea 
PA a 2 ee eee 
ht Es a EO st Ui 1 sp EG et a a UV Ds st os Dh eel Ui hs OP cs ce 

{ 

for (m=(i+sf_size); m<(i+2*sf_ size); m++) 


{ 


l=1+sf_size+l1; 


1f((sf_size==2) && (fF L{[1]!'="*b") & & (2 Ai) bate) & & 
(f_1[l]!='d’) es (f_l[l]!='‘e") & & (Ey Llane") & & 
(22) [aopeha Go) see CE a Ed J let) 

tlm) ="a’ 3 | 
else if ((sf_size==4) &&(f_l[{l]!s’c’) &&(f_l[l]!='d’) &«& 
¢oe a pee ce (£1 bp tate) &&(f_l[l]!='g’) & & 
coe (es es oe 

see ee 0d ise ope, 


else if ((sf_size==8) && (f_l[l]!='d') && ( 
(EL i): ee (frie gg") 68 (F211) fret 
Fhimps ers 
else if ((sf_size==16) && (f_l[l]!=#‘e’) s&s (E-file e") 2G 
(fF_l[l]!='’g’) && (f_1[1]!="h')) 
fl[mj=‘d'; 
else if ((sf_size==32) && (f_l[l)!='f') && ET tee) ee 
oa ea a) eee CE ara) ee oe aa aay 
f_l[{m]j=’e’; 
else if ((sf_size==64) && (eel ero). eee: ea ere ae 
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(rea ee oe | et ae 
££ i[m]=’f’; 
else if ((sf_size==128) && (f_lf{lJ!='h’) && (f_1[1l]!='i') 
aa Saamme Os a 
f1l[{mj=’g'; | 
else if ((sf_size==256)&&(f_l[1l]!='i')&(f_1[1l]!='j")) 
f l[{m]j=‘h’; 
else if ((sf_size==512)&&(f_lo[1] !='j’)) 
f_i[m])=’i’; 
else if (sf_size==1024) 
feiss 
} 
} 
} 


1=1+2*sf size; 


} 


f_1l[pow(2,num_var)]=’\0’; 
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/*- Prepare result to be converted into binary tree ----- xf 


LOr (12=0>7 1=(pow(2,num var)=1)> i++) 
f_vector[iJ='*’; 
for (1=(pow(2,num_var)-1); i<(pow(2, (num_var+1))+1); i+4) 


{ 


Lf (1==(powt2, (num vare1)))) 
{ 

f vector[i]J='\0’; 

break; 
} 


else 


1f (£_1[1-(pow(2,num_var)-1)]=='0') /* Lowest level 
subfunction */ 
£-Vector|[1. |e" 0"? 
if (f_1[1i-(pow(2,num_var)-1)]=='1') 
E-Vecror (3: |=*L" 
1f (f_1[i-(pow(2,num_var)-1)]==’a’) 
f vector[i]="'a’"; 
if (f_l1[i-(pow(2,num_var)-1)]==’b’) 
f vector[i]=’b’; 
1f (f_1[i-(pow(2,num_var)-1)]==’c’) 
t -vector[i}="c" > 
1f (f_1[i-(pow(2,num_var) -1) ]==’d’) 
f vector[i]J=‘d’; 


1f (£_1[1-(pow(2,num_var)-1)]==’e’) 
f vector[ij=’e’; 

if (f£_1[1i-(pow(2,num_var) -1)]==’'f') 
b vector[i)/="f"% 


if (£_1[1-(pow(2,num_var)-1)]=='g 
fo VeOCrcor (2 |.=' qs 

1f (f_1[1i-(pow(2,num_var) -1)]=="'h’) 
f vector[ijJ=’h’'; 

1f (f_l1[i-(pow(2,num_var)-1)]=='i’) 
f_vector[ijJ=’i'; 

if (f_1[1i-(pow(2,num_var)-1)]=='j’) 

/* Highest level subfunction */ 
£F vector[iJ='3’; 


} 
jf PES = Create Binary Tree-------~---------------------- a 
root=create_tree(f_vector,0, (pow(2, (num_var+1) )-1) ); 


for (1<0+ a<num: vare 144) 
preorder (root, &node_remain) ; /* Call. -Lunction. to 
Simplify binary tree 
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a aa OUEDUG eh eetasoeeeseaseeSsoee i] 


cross _flag=0; 

inorder(root, &cross_flag) ; /* print out all nodes in 
order & do crossing check */ 

node_remain=2; 


*node_cnt=2; 


/* Count the number of nodes & free memory used by binary 
tree-- */ 


postorder(p, root, node_cnt, num_var); 


1f (cross_flag==2) 
{ 
printf£("\n Planar BDD\n "); 
for (j=0; j<num_var; j++) 
Eprintf(ofp," %d",x_order [num _var-j-1]) + 


Eprintft (ofp, * td d\n “,threshold, *node_cnt) ; 
if (min_node>*node_cnt) 
min_node=*node_cnt; /* seek minimum nodes BDD sa 6 
print tree(p, root ,num_var, x_order); 
} 
if (cross_flag==4) 

{ 

if (min_node>*node_cnt) 

min_node=*node_cnt; /* seek minimum nodes BDD */ 

printf£("\nCrossing occurs "); 


[taser = Write to file if required---------~---------.-.- st A 


/* Ascending order */ 
1f (num_var==4) 
{ 
if ( (x_order[0]>x_order [1] ) && (x_order [1]>x_order[2]) 
&& (x_order [2]>x_order[3])) 
print_to_file(*node_cnt, p, root, num_var, filename, 
threshold,w_in, x_order) ; 


} 


if (num_var==5) 


{ 
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a 


1 
((x_order[0]>x_order[1]) &&(x spears order [2])&&(x_order [ 
2]>x_order[3]) &&(x_order[3]>x_order[4]) 

print_ to_file(*node_cnt, p, ee num_var, filename, 
threshold,w_in, x_order) ; 


} 


if (num_var==6) 

{ 

1 f 
((x_order[0]>x_order([1]) &&(x_order[1]>x_ sgpeagedar && (x_order [ 
2])>x_order[3])&&(x_order[{3]>x_order[4]) &&(x _order [4] >x_order 


ie) 
print_ to_file(*node_cnt, p, root, num_var, filename, 


threshold,w_in, x_order); 


} 


if (num_var==7) 

‘ 

a f 
((x_order[0]>x_order [1] ) &&(x_order[1]>x_order[2]) &&(x_order [| 
2]>x_order[3])&&(x_order[3]>x_order[4]) )&& (x_order [4]>x_order 


[5] ) &&(x_order [5]>x_order[6]) ) 
print_ to_file(*node_cnt, p, root, num_var, filename, 


threshold,w_in, x_order) ; 


} 


/* Descending order */ 


1f (num_var==4) 


{ 


si ions Paella a 
r[2]<x_order [3] ) && (x_order [3]>x_order[4]) 

print_ to_file(*node_cnt, p, eel num_var, filename, 
threshold,w_in, x_order); 


} 


1f (num_var=<=5) 


{ 


oy Epes basepair Saved tegen i i ata cacy 
r{[2]<x_order[3])&&(x _order [3]<x_order[4]) 

print_ to_file(*node_cnt, p, het num_var, filename, 
threshold,w_in, x_order); 
} 


if (num_var==6) 


{ 


68 








1f ( (x_order [0]<x_order[1]) &&(x_order [1]<x_order[2]) &&(x_orde 
r[2]<x_order[3])&&(x_order [3]<x_order [4] ) &&(x_order[4]<x_ord 
er i))-) 

print_to_file(*node_cnt, p, root, num_var, filename, 
threshold,w_in, x_order); 


} 


1f (num_var==7) 

1 si 
((x_order [0]<x_order [1] ) && (x_order [1] <x_order[2])&&(x_order [ 
2]<x_order [3] ) && (x_order [3] <x_order [4] ) & (x_order [4]<x_order 
[5] ) && (x_order [5]<x_order[6])) 

print.to.file(*node cnt; p, root, num var; filename, 
threshold,w_in, x_order) ; 


} 


} /* if a permutation is found */ 


} while (carry[num_var]==0); /* while permutation 
is not completed */ 
fprintf(ofp, "\n Minimal BDD has $d nodes \n",min_node) ; 
tclose(ofp) ; 
} /* fibonacci Loop , from 1 to T_max */ 


/* fclose(ofp); */ 


} /* end of main */ 


/*create a linked binary tree from an array */ 
[EEE ERICA SAT RI RK ARR KR Ae i Se Rte ie Be He Kea ek te ae] 
BTREE create_tree(DATA a[], int i, int size) 
{ 
if (i >= size) 
return NULL; 
else 
return(init_node(afi], i, create_tree(a, 2*i+1, size), 
create_tree(a, 2*1+2, size))); 


} 
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[BREE RER REE EE EAR ALI AR HA RY 


/* Creating a binary tree */ 
f REALE RR RAR ARIE R BA Ae eK 


BTREE new_node() 

{ 

return (malloc(sizeof(NODE))); 
} 


BTREE init_node (DATA dl, int i, BTREE pl, BTREE p2) 
{ 

BREE. 

t=new_node(); 

t->index=i+1; 

t->d=dl1; 

t=slertap.. 

tEaSragntap2: 

LEC: "e 


void inorder(BTREE root, int *cross_flag_ptr) 
{ 

1f (root !=NULL) 

{ 

inorder (root->left, cross_flag_ptr) ; 

cross_test (root->d,cross_flag_ptr) ; 

/* prance’ (“$e",. root <sd)s */ 

inorder (root->right,cross_flag_ptr); 


[ Re Scan for mergeable nodes --------~-------------_- aa 
void preorder (BTREE root, int *node_remain) 


1£ (root !=NULL) 
{ 
merge(root, &node_remain) ; 
preorder (root->left, &node_remain) ; 
preorder (root->right, &node_remain) ; 
} 

} 


/*-~ count the number of nodes left in simplified BDD ----*/ 


void postorder(DATA p[], BTREE root, int *node_cnt, int 
num_var) 


{ 
DATA rd; 
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discount=0; 


oy 


Tit 
1£ (root !=NULL) 


{ 


postorder(&p[0], 
postorder(&p[0], 


root->left, 
root->right, 


p[root->index] =root->d:; 


rd=root->right->d; 
1=root->right->index; 


1£f ((root->d)=='*") 








*node_cnt = *node_cnt+l1 ; 


L£ 
{ 





node-cnt., 
node_cnt, 


((rd!='0')&&(rd!='1')&&(rd!i="*"')) 





num_var) ; 
num_var) ; 


1f( (root->left->d=='0') | | (root->left->d=='1')) 


discount=1; 


1f£((i>(pow(2, (num_var-1) ) -1) && (i<pow(2,num_var) )&&(rd=='a’))) 


discount=1; 


1f((1>(pow(2, (num_var-2))-1)&& (i<pow(2,num_var-1) ) &&(rd=='b’ 


4) 
discounte=1; 


1f((1>(pow(2, (num_var-3) )-1) && (i<pow(2,num_var-2))&&(rd=='c! 


Pp 
discount=1; 


1f((1>(pow(2, (num_var-4) )-1) && (i<pow(2,num_var-3) ) && (rd== 


ry 
discounte=1; 


al 
Pe 
discount=1; 


((1>(pow(2, (num_var-5) )-1) &&(i<pow(2,num_var-4) )&&(rd==-’e! 


1f((1>(pow(2, (num_var-6) )-1) &&(i<pow(2,num_var-5) )&&(rd=='£! 


Pe) 


discount=1; 


1f((1>(pow(2, (num_var-7) )-1) && (i<pow(2,num_var-6) ) && (rd=="g' 


))) 
discount=1; 


1f((1>(pow(2, (num_var-8) ) -1) && (i<pow(2,num_var-7))&&(rd=='h’ 
)) 


discount=1; 


1f((1>(pow(2, (num_var-9) ) -1) && (i<pow(2,num_var-8))&&(rd=='4i! 


pee 


))) 
discount=1; 


1f((1>(pow(2, (num_var-10) )-1) && (i<pow(2,num_var-9)) &&(rd=='3 
rays) | 
discount=1; 
1f (discount==1) 


{ 


*node_cnt=*node_cnt-1; 
discount=0; 
} 

} 


Free(root) ; 


} 


[TEES Tees SMO Lt eat ton LOuLInNeS Tor Binary Tree *<*> */ 


/* Return a ‘0/1’ if left branch & right branch are both 
‘0/1’ or else mark node with ‘*’ to indicate that merging is 


not possible */ 


DATA merge (BTREE node_ptr, int *node_remain_ptr) 
{ 

DATA t,1n,rn; 
ln=node_ptr->left->d; 
rn=node_ptr->right->d; 


if ((in==rn) && (node_ptr->index! =1) ) 
{ 
Lt 
({in=='0')// (In=='1') 1] (In==’a’) || (1n=='b’)/[ (1In=='c"’) || (1n= 
=‘d‘)) 
{ 
t=node_ptr->left->d; 


node_ptr->d=t; 
node_ptr->left->d=' '’; /* terminate branch */ 


node_ptr->right->d=’ ‘'; 

return. E; 

} 

else if ({in==’e’) || (In==’f') |! (ln=='’g’) || (lns=’h’) 
ae me C2 aa a 

{ 


t=node_ptr->left->d; 


node_ptr->d=t; 
node_ptr->left->d=' ‘; /* terminate branch */ 
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/ 


node_ptr->right->d='’ '; 
Secure 


[Reese Sa eee Print result ------------------- 


void print_tree(DATA p[], BTREE root, int num_var, 


x_order[]) 


Ine 245 

DrIneet" an"); 
printf ("X%d 

eC" 7% order [0)],p [i] )% 


prince’ ("\n\n*)-; 


Prancr ("xsd oc", x order [1], 0 ([2]); 
Priner ( SC "7p lol); 
printf ( 


sara a 6 Ms ) : 
DrINtCL ("xsd sc", x order [2],p[4]); 
for (9257 5J<8. ja) 
PYinte (* SC, DIL) 


1£ (num_var>2) 


{ 


Drier nk )¢ 
1f (num_var!=3) 
printf ("xd ",%order([3}); 
printt ("sa",p[8)]): 
For (3=9; 3<16; 3++) 
Perntet™ Se" >: “PAI 4 


if (num_var>3) 
{ 
Pprianctt("\n\n") 
1f (num_var!=4) 
Drintt(*xXtd “yx order [4] ); 


printf ("%c",p[16]); 
for (9=172 3<323 544) 
brintr(” sc",pi(jl]): 
if (num_var>4) 


{ 
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int 








Prine a 3 )% 
Lor (j=322% 3264. 5 ae) 
Print 26", 015 173 


RS Si eee Print result to file---------------~-------- at 


void print_to_file(int node_cnt, DATA p[{], BTREE root, int 


num var, 
char filename[], int threshold, int w_in[], int x_order[]) 


{ 
Ine: Ae 
FILE *ofp; 


ofp = fopen("D_perm.7vd","a"); 


fprintf(ofp,"\nThreshold function F=("); 
bor ( 220% 22nuit. vars 464) 
ERrIMNtE (OL p. "Sd." we in 7] ) 
fprintf(ofp,"; dad ) \n",threshold); 


forintt(ofp,"Order evs 

bor | 1203 a<ntnivar* 744) 

Eprantt (ofp, "xed. “x order li} \- 
fprintf(ofp,"\nNumber of nodes = %d \n\n",node_cnt) ; 


[ERKKKKKKKKKKKRKKKKKK Lureral SF a EA We Tee is eee De sR Ie I EEA oR ade 


fprintf(ofp, “\Binary Tree Map\n"); 
for (1=64; i<(pow(2, (num_var+1))) ; i++) 


{ 
Eee 
fprintf( ofp, "Index = $d SC. Ase pla 
} 
Pprincr (ote... “\nin')> 


[RR EKKKKKKKKEKKK Literal DUNE Re Fete ieee eR tee Me cae Nee OR NR ee aR ae 


[® moma e --- ~~ level 1 */ 
Pprintt (orp; "\nin") = 
EPrianer (loco <" 

SC", Dit): 
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Eprinte (ofp. Wn in) 
PPrintt (ofo." 
LDEIAncrEe orp, « 


oo 


CPO |). 


1f (num_var>1) 

{ 

EDYINCE OLD; Vii") 3 

fprintf (ofp, " oc", p[4]); 

FOr (9=5; 97288 “jee 

forint’ (orp, 5 oleae © be i ey 
} 


1f (num_var==2) 
Pporintt (ofp, "an" )-3 


[* mmo n-ne ~~ ~~ =~ - level 4 */ 


1f (num_var>2) 
{ 
Fprintt (orp, n\n") 
/* if (num_var!=4) 
EPrintr (orp, "xs "7x Order [4] ) a 
Eprintt(ofo," sc ", p[8]); 
FOr (729% 9eLey: 34) 
eprint (ofp, = oC" Dial) 
} 


1f (num_var==3) 
EPLinet (ofp, win") 


[Ro mmm mn nnn nnn == ~~ - Level. 5. %/ 
if (num_var>3) 

{ 

EDrinktt (ofpy Vivi) 

Pprintet (orp." $c", p[16]); 


POE (jel ya32s jae) 
fprintt (ofp;* SC°4 (Did) 3 


} 


if (num_var=<=4) 
fprintt (ofp, “\n\n"):: 


[RB mom mw nnn nnn nn ne level 6 *Y/ 


1f (num_var>4) 


fis 





{ 


Porinte (Op, \n in): 
fprinet (ofp. Be", 137) \s 


£Oe (jas 37 7264" 744) 
COrInet (Orp,” Louden oy aa 


} 


1f (num_var==5) 
Porn L (Ole. nin 


[*® --------~~--~------------------ level 7 */ 


1f (num_var>5) 
{ : 
Pprince (ofp, nin") 

Lor (j=64> 3<128? 544) 
EOriIntt (6fp, "se "3015144 


Epraintt (oto, \n\n\n\n" )4 
} 


Eclose(ofp); 
} 
[*® ~--------- BxXaming if there 1s Crossing===—-s=42-4-== a 4 
cross_flag =0 : no ‘0’ or ‘1’ occur yet 
cross_flag =1 : first ‘0’ occured 
cross_flag =2 : first ‘1’ occured 
cross_flag >=4 : >1 ‘'0->1’ transition has occured 
Tr a wr a se eee ee ee * / 
void cross_test (DATA d, int *cross_flag) 
{ 
1f ((*cross_flag==0) &&(d=='*")) 
*cross_flag=0; 
1f ((*cross_flag==0) &&(d=='0’)) 
*cross_flag=1; 
1f ((*cross_flag==1)&&(d=='1")) /* Farst crossing occurs 
all 


the time */ 
*cross_flag=2; 


if ((*cross_flag==2) &&(d=='0')) 
*CLOsst1ags3 > 
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1£ ((*cross_flag==3) &&(d=='1')) y= => 1 Crossing occurs 
ae 


} 


*cross.flag=4; 


[RRR KKK KKK KKK KKK KEK KK KK RK KEK KK KK KEK KK RK KEK KEK KKKKKKEKKKKKKKK KKK 


Compare sub-function to see that sfl is not all '0’ or all ‘1’ 
before EOF; return an integer 0 if the above is true, else 
return an integer 1 return -1 is error calling of function 


REE AS Se a re Me RR aaa Re ee IIE eR ee Ke EDS Ie eee ie Dn dee ae ea SRY 


int sfl_cmp(DATA sf1[]) 
{ 


int num_char=0,j=0; 
int ent01=0; 


while (sf1[j+1]!=EOF) 


i1f((sf1[j]==sf1[j3+1]) &&(sf1[j+1] !=EOF) &&((sfl[j]=='0') 1] (sf1 
[j]=='1"))) 
cnt01l++; 
seee 


} 
while (sf1l[num_char] !=EOF) 
num_char++; 
return (num_char-cnt01); 


} /* end of function */ 


fa 
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